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ABSTRACT 

An accurate knowledge of the mineralogy (chemical composition and crystal structure) 
of the silicate dust in the interstellar medium (ISM) is crucial for understanding its 
origin in evolved stars, the physical and chemical processing in the ISM, and its sub- 
sequent incorporation into protostellar nebulae, protoplanetary disks and cometary 
nuclei where it is subjected to further processing. While an appreciable fraction of 
silicate dust in evolved stars, in protoplanetary disks around pre-main sequence stars, 
in debris disks around main sequence stars, and in cometary nuclei is found to be 
in crystalline form, very recent infrared spectroscopic studies of the dust along the 
sightline toward the Galactic center source Sgr A* placed an upper limit of ~ 1.1% on 
the silicate crystalline fraction, well below the previous estimates of ~ 5% or ~ 60% 
derived from the observed 10 /mi absorption profile for the local ISM toward Cyg OB2 
No. 12. Since the sightline toward Sgr A* contains molecular cloud materials as revealed 
by the detection of the 3.1 and 6.0 /im water ice absorption features, we argue that by 
taking into account the presence of ice mantles on silicate cores, the upper limit on 
the degree of silicate crystallinity in the ISM is increased to ~ 3-5%. 

Key words: ISM: dust, extinction - Galaxy: center - infrared: ISM: line and bands 
- ISM: diffuse ISM, molecular clouds 



1 INTRODUCTION 

The mineralogical composition of dust contains important 
information about its origin and evolution, and may reveal 
the physical, chemical and evolutionary properties of the 
astrophysical regions where the dust is found. Silicate, one 
of the major components of cosmic dust species, is ubiqui- 
tously seen in various astrophysical environments, ranging 
from the Galactic diffuse interstellar medium (ISM) , HII re- 
gions, and the dust torus around active galactic nuclei to 
dust envelopes around evolved stars, protoplanetary disks 
around pre-main sequence stars, debris disks around main 
sequence stars, cometary comae and interplanetary space. 
Their chemical composition and crystal structure vary with 
local environments. 

Interstellar spectroscopy provides the most diagnos- 
tic information on dust composition. In the infrared (IR), 
the strongest interstellar spectral features are the 9.7 and 
18 fim absorption (or emission) bands which are generally at- 
tributed to the Si-O stretching and O-Si-0 bending modes, 
respectively. While the absorption profiles measured in labo- 
ratory for crystalline olivine and pyroxene show many sharp 
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features, in the diffuse ISM, the observed 9.7 and 18 /im sil- 
icate features are broad and relatively featureless, suggest- 
ing that interstellar silicates are mainly amorphous rather 
than crystalline. Since the crystallinity of silicate dust is inti- 
mately connected to the energetic processing occurred in the 
evolutionary life cycle of dust and depends on the environ- 
mental properties, of particular interest to silicate astromin- 
eralogy is the abundance of crystalline silicates in different 
astrophysical environments, especially in the diffuse ISM. 

Recently, a number of studies have been made to de- 
termine the crystallinity degree of silicates in the diffuse 
ISM. From the observed 9.7 fim absorption profile for dust 
in the local ISM toward Cyg OB2 No. 12, Li & Draine (2001) 
estimated the fraction of Si in crystalline silicates to be 
<5%. Demyk et al. (1999) derived an upper limit of ~ 1— 
2% of crystalline silicates in mass toward two massive pro- 
tostars. However, Bowey & Adamson (2002) argued that a 
complex mixture of crystalline silicates (60% by mass) and 
amorphous silicates (40% by mass) could explain the ob- 
served smooth silicate absorption profile at 9.7 /im toward 
Cyg OB2 No. 12. More recently, Kemper, Vriend & Tielens 
(2004) placed a more strict constraint of < 1.1% on the crys- 
tallinity of interstellar silicates, based on a detailed analysis 
of the 9.7 /im feature obtained with the Short Wavelength 



2 M. P. Li, G. Zhao & A. Li 



Spectrometer (SWS) on board the Infrared Space Observa- 
tory (ISO). This upper limit ~ 1.1% of crystalline silicates 
was derived from a direct comparison of the Sgr A* spectrum 
with theoretical spectra for pure silicates. 

However, it is evident from the detection of the 3.1 and 
6.0 fim water ice features respectively attributed to the O-H 
stretching and bending modes that there are molecular cloud 
materials along the line of sight toward Sgr A* (McFadzean 
et al. 1989; Tielens et al. 1996; Chiar et al. 2000). This is 
further confirmed by the detection of solid CO2 absorption 
toward Sgr A* (Lutz et al. 1996; de Graauw et al. 1996). 
In order to infer the precise mineralogical composition of 
dust, we have to take the grain ice mantles into account 
when interpreting the observed silicate features. Indeed, in 
this paper we show that ignoring the ice mantles coated on 
the silicate cores can result in an underestimation of the 
crystalline degree of silicate dust, as much as ~ 3-5%. 

The purpose of this work is to address the question 
how the ice mantles affect the determination of dust min- 
eralogical composition (i.e. whether and to what degree the 
inclusion of ice mantles will hide the sharp features of crys- 
talline silicates and lead to an underestimation of the sil- 
icate crystallinity) , not to model any specific astronomical 
objects. Therefore, neither the specific choice of silicate di- 
electric functions nor the precise constituents of the ice man- 
tles would affect our conclusion. 
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Figure 1. Absorption cross sections of spherical grains for 
Vi ce /V s ii = 0.55. Dotted: amorphous silicate core-ice mantle 
grains (CJJJ); dot-dashed: crystalline olivine silicate core-ice man- 
tle grains {C^, B ); long-dashed: equal- volume ice spheres for the 
ice mantles; "solid: C abs = + AYC^ where AY = 5% 

is the crystalline olivine silicate fraction; short-dashed: C a b s = 
C™ + AYC^ S with AY = 3%. Note that water ice produces 
a shoulder at A > llitm. This should not be interpreted as a 
crystalline silicate feature. 



2 VOLUME FRACTION OF WATER ICE 

The sightline toward the Galactic center source Sgr A* suf- 
fers about ~ 30 mag of visual extinction (e.g. see McFadzean 
et al. 1989), to which molecular clouds may contribute as 
much as ~ 10 mag (Whittet et al. 1997). Except for H 2 
features, there are many additional molecular absorption 
features, attributed to C0 2 , NH 3 , CO, CH 3 OH, CH 4 and 
other species. These icy molecules might accrete on the pre- 
existing silicate core inside dense molecular clouds and form 
an ice mantle. As for further evolution, the accreted icy grain 
mantles could be photoprocessed and converted into organic 
refractory residues (Greenberg et al. 1995). 

Assuming the 3.1 /jm water ice absorption results from 
the water ice mantles uniformly coated on silicate grains 
which produce the 9.7/xm absorption, we estimate the vol- 
ume ratio of the ice mantles to the silicate cores from the 
observed optical depths - T3.1 for the 3.1 /im O-H feature, 
and T9.7 for the 9.7 /^m Si-O feature. Let V B n and Vice re- 
spectively be the volumes of the silicate cores and the ice 
mantles. Let Af c l and be the integrated band strength 
for the 3.1 /im O-H feature and the 9.7 /im Si-0 feature, 
respectively. For spherical silicate core-ice mantle grains, we 
use Mie theory (Bohren & Huffman 1983) to obtain Af c * = 
Lx„Jf%£lV)dX and AH = J 9 , ril JCH/V)dX, where 
Cice (C'sii) i s the continuum-subtracted absorption cross 
section of the 3.1 (9.7 (im) ice (silicate) feature, and V is 
the volume of a spherical grain with radius a. For the Sgr A* 
sightline, r 3 .i « 0.50 and r 9 . 7 « 3.52 (Chiar et al. 2000). We 
therefore obtain Vicc/Vk = (t 3 .i/t 9 . 7 ) (A%i /Af£) as 0.55, 



assuming that all silicate grains are coated by a layer of ice 
mantle rj 



3 SPHERICAL GRAINS 

We consider 3 types of dust materials: amorphous olivine, 
crystalline olivine, and pure water ice. We adopt the di- 
electric functions of Dorschner et al. (1995) for amorphous 
olivine, of Li & Draine (2001) for crystalline olivine, and of 
Hudgins et al. (1993) for water ice. 

We first consider spherical silicate core- water ice mantle 
grains. We take the silicate core size to be a s n = 0.1 fj,m, a 
typical size for interstellar dust. We should note that the 
choice of an exact grain size is not critical since in the 
wavelength range considered here submicron-sized interstel- 
lar grains are in the Rayleigh regime. Also because of this, 
we do not need to consider dust size distributions. 

In Figure [T] we plot the absorption cross sections for 
amorphous olivine silicate core-H2 ice mantle grains and 
for crystalline olivine silicate core-H20 ice mantle grains. 
As shown in Figure [T] the absence of narrow features near 
10.0 and 11.1 /zm in the absorption spectrum suggests that 

In reality, a thicker ice mantle would be expected for the silicate 
dust in the molecular cloud component which contributes about 
10 mag to the total ~ 30 mag of visual extinction toward Sgr A* . 
An ice mantle as thick as Vi ce /V a n ~ 7.6 was implied for the dust 
in the dense molecular cloud toward the field star Taurus Elias 16 
(Bowey et al. 1998). A thicker ice mantle would hide up a larger 
fraction of crystalline silicates. 
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the inclusion of water ice mantles would hide up ~3-5% 
crystalline mass fraction without being noticed. Note that 
water ice has a broad absorption band at ~ 12.2 /im (see 
Fig.[[] Leger et al. 1983; Hudgins et al. 1993). The inclusion 
of an ice mantle on silicate dust (no matter it is amorphous 
silicate or crystalline silicate) results in a weak shoulder at 
A > 11 fim. This is noticeable even in the model spectrum 
of ice-coated pure amorphous silicate dust (see Fig.0. One 
should caution that this should not be interpreted as a crys- 
talline silicate feature. 

We also consider a model in which only the molecular 
cloud dust is coated with a layer of ice mantle. If we assume 
that silicate dust and carbon dust equally contribute to the 
~ 10 mag visual extinction, the volume ratio of the ice man- 
tle to the silicate core for the molecular cloud dust would be 
Vice/Kiii ~ 0.83. The bare diffuse cloud dust is responsible 
for the remaining ~ 20 mag visual extinction. In Figure [2] we 
show the model cross sections obtained by summing up 1/3 
of that produced by the molecular cloud dust and 2/3 of 
that produced by the diffuse cloud dust. Similarly, ~3-5% 
crystalline silicate would be hidden because of the presence 
of an ice mantle on the molecular cloud dust toward Sgr A* Q 

So far, we have only considered olivine silicate dust. It 
is possible that there is a considerable amount of pyroxene 
dust in the ISM (see Bowey & Adamson 2002 and refer- 
ences therein). For illustration, we carry out similar calcula- 
tions for pyroxene and obtain similar conclusions (see Fig. [2]). 
Again, we should stress that the purpose of this work is not 
to model any specific astronomical objects, but to investi- 
gate the effects of ice mantles on the silicate crystallinity es- 
timation. We do not expect either olivine or pyroxene alone 
could closely reproduce the interstellar silicate absorption 
feature observed toward Sgr A*. For comparison, we also 
show in Figure[2]the ISO spectra of Sgr A* obtained by Kem- 
per et al. (2004) and Gibb et al. (2OO4)0The fact that, while 
the 10 /jm feature of olivine peaks at relatively longer wave- 
lengths than the observed spectra, pyroxene peaks at too 
short wavelengths, suggests the co-existence of both olivine 
and pyroxene in the ISM. 



4 SPHEROIDAL GRAINS 

However, interstellar grains must be nonspherical as in- 
dicated by interstellar polarization. We now consider 



2 Admittedly, the approach adopted here is simplified. The 
molecular cloud silicate dust may differ from that of the diffuse 
cloud (e.g. their 10 fim Si— O features may have different line pro- 
files [particularly line widths], see Bowey, Adamson, & Whittet 
2001 and references therein). Their optical properties which are 
temperature-dependent (Bowey et al. 2001) may also differ from 
each other since the molecular cloud dust is generally colder than 
the diffuse cloud dust (e.g. see Greenberg & Li 1996a). Moreover, 
it is not clear if the molecular cloud silicate dust has the same 
iron fraction as the diffuse cloud silicate dust (it is well recognized 
that, as most recently experimentally demonstrated by Bowey et 
al. [2007], the optical properties of silicate dust are very sensitive 
to the proportion of iron). 

3 We note that the ISO spectrum of Sgr A* appears to have a 
feature at ~6.9/xm, attributed to crystalline melilite (Bowey & 
Hofmeister 2005). 
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Figure 2. Comparison of the ISO spectra of Sgr A* (thick solid: 
Gibb et al. 2004; thick dashed: Kemper et al. 2004) and the ab- 
sorption spectra of olivine (thin solid) and pyroxene (dot-long 
dashed) spheres with 3% crystallinity, with 1/3 of the dust coated 
by a layer of ice mantle (Vi co /V^ii = 0.83) and 2/3 of the dust be- 
ing bare silicates. The A > 11 fim shoulder of the model spectra 
arises from the ice mantles (also see Fig. [I}. Also shown is the 
absorption spectrum of pure amorphous olivine dust without an 
ice mantle (dot-short dashed). 



spheroidal shape which can be exactly solved in the Rayleigh 
limit (see Li et al. 2002 and references therein). 

We consider confocal silicate core-ice mantle spheroidal 
grains with a volume-equivalent sphere radius r^J} = 

(ab 2 ) 1 ^ 3 = 0.1 ^m, where a and 6 are the semi-axis along and 
perpendicular to the symmetry axis, respectively. Again, we 
take the volume ratio of the ice mantles to the silicate cores 
to be Vico/Vsii = 0.55. 

We first consider prolate grains with r^J = 0.1 /im for 
the silicate cores and a/b = 3 for the ice mantles0 In Fig- 
ure^ we plot the absorption cross sections for ice coated- 
amorphous silicate prolates and ice coated-crystalline pro- 
lates. Just like spherical grains, the inclusion of water ice 
mantles would hide up ~ 3% crystalline mass fraction with- 
out being noticed. Similarly, we perform the same calcula- 
tions but for oblate grains with a/b — 1/2 for ice mantles[f| 
As seen in Figure [3]3, a crystallinity degree of as much as 
~ 3% can be covered by the ice mantles. 



4 Greenberg & Li (1996b) found that prolates of a/b = 3 provide 
an almost prefect match to the 10, 18 (im silicate polarization 
features of the Becklin-Ncugebauer (BN) object. 

5 Lee & Draine (1985) and Hildebrand & Dragovan (1995) found 
that a/b = 1/2 oblates fit well the 3.1 fim ice polarization and the 
10 fim silicate polarization of the BN object. 
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Figure 3. Same as Figure[T]but for (a) a/b = 3 prolates and (b) 
a/b = 1/2 oblates. 



5 GRAINS WITH DISTRIBUTIONS OF 
ELLIPSOIDAL SHAPES 

The shape of a spheroidal grain is characterized by its eccen- 
tricity e (and the depolarization factors L", L x via eq. [4] of 
Li et al. 2002) . In §4 we consider grains with a single shape. 
We now consider an average of a distribution of shapes. We 
assume two kinds of shape distribution functions: 1) dP/dL^ 
= constant, i.e., all shapes are equally probable (Bohren & 
Huffman 1983); 2) dP/dL 11 = 12L 11 [1 - L 11 ] 2 (Ossenkopf, 
Henning, & Mathis 1992)0 

Averaging over the shape distribution, we 
have the resultant absorption cross section C a b s = 
J dL"dP/dL"C a bs(L"), where C a bs(L") is the absorption 
cross section of a particular shape L" (note L" is for the 
mantle; the core depolarization factor is derived from 
eqs.[4-6] of Li et al. 2002). Again, we assume confocal 
geometry for core-mantle grains, with the above dP/dX" 
(as well as e) applying to the outer surface. 

The results are shown in Figure [4] Similar to spherical 
grains and grains of a single ellipsoidal shape, with ~ 3% 
crystalline silicates included, the absorption profiles still do 
not seem to exhibit the sharp features of crystalline silicates. 



6 SUMMARY 

We have investigated the effects of ice mantles coated on sili- 
cate cores of various shapes on the determination of the crys- 

6 This distribution peaks at spheres (LH = L ± =l/3) and is sym- 
metric about spheres with respect to eccentricity e. It drops to 
zero for the extreme cases: for infinitely thin needles (a ^> b) 
L" — > 0; for infinitely flattened pancakes (a <C b) L" — * 1. 



Figure 4. Same as Figure [T] but (a) for grains with a uniform 
distribution of ellipsoidal shapes (i.e. dP/dL " =constant) and (b) 
for grains with dP/dlJ = 12.ZJ [1 - L^] 2 . 

tallinity degree of silicates. It is found that for the dust in the 
line of sight toward the Galactic center source Sgr A* , ~ 3% 
crystalline silicates could be hidden by ice mantles, well ex- 
ceeding the upper limit ~ 1.1% of Kemper et al. (2004) de- 
rived from the assumption of the absence of dense molecular 
materials in this line of sight. We take the standard approach 
by assuming that the silicates are amorphous and then 
adding in crystalline components. But as shown in Kemper 
et al. (2004), the smooth, featureless 10 /im amorphous sili- 
cate spectrum allows as much as ~ 2.2% crystalline silicates 
(or even higher; see Bowey & Adamson 2002). Therefore, 
the total allowable degree of crystallinity would be ~ 5%, 
consistent with the earlier estimates of Li & Draine (2001). 
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